Dual Quantum Monte Carlo Algorithm for Hardcore Bosons

We derive the exact dual representation of the bosonic Hubbard model which takes the form of conserved current loops. The hardcore limit, which corresponds to the quantum spin-${1\over 2}$ Heisenberg antiferromagnet, is also obtained. In this limit, the dual partition function takes a particularly simple form which is very amenable to numerical simulations. In addition to the usual quantities that we can measure (energy, density-density correlation function and superfluid density) we can with this new algorithm measure efficiently the order parameter correlation function, $, |i-j|\ge 1$. We demonstrate this with numerical tests in one dimension.Comment: 15 pages, 4 figures . Talk given at CCP1998, Granada, Spai

Langevin Simulations of a Long Range Electron Phonon Model

We present a Quantum Monte Carlo (QMC) study, based on the Langevin equation, of a Hamiltonian describing electrons coupled to phonon degrees of freedom. The bosonic part of the action helps control the variation of the field in imaginary time. As a consequence, the iterative conjugate gradient solution of the fermionic action, which depends on the boson coordinates, converges more rapidly than in the case of electron-electron interactions, such as the Hubbard Hamiltonian. Fourier Acceleration is shown to be a crucial ingredient in reducing the equilibration and autocorrelation times. After describing and benchmarking the method, we present results for the phase diagram focusing on the range of the electron-phonon interaction. We delineate the regions of charge density wave formation from those in which the fermion density is inhomogeneous, caused by phase separation. We show that the Langevin approach is more efficient than the Determinant QMC method for lattice sizes $N \gtrsim 8 \times 8$ and that it therefore opens a potential path to problems including, for example, charge order in the 3D Holstein model

Exact duality and dual Monte-Carlo simulation for the Bosonic Hubbard model

We derive the exact dual to the Bosonic Hubbard model. The dual variables take the form of conserved current loops (local and global). Previously this has been done only for the very soft core model at very high density. No such approximations are made here. In particular, the dual of the hard core model is shown to have a very simple form which is then used to construct an efficient Monte Carlo algorithm which is quite similar to the World Line algorithm but with some important differences. For example, with this algorithm we can measure easily the correlation function of the order parameter (Green function), a quantity which is extremely difficult to measure with the standard World Line algorithm. We demonstrate the algorithm for the one and two dimensional hardcore Bosonic Hubbard models. We present new results especially for the Green function and zero mode filling fraction in the two dimensional hardcore model.Comment: 14 pages, 15 figures include

Haldane phase in the sawtooth lattice: Edge states, entanglement spectrum and the flat band

Using density matrix renormalization group numerical calculations, we study the phase diagram of the half filled Bose-Hubbard system in the sawtooth lattice with strong frustration in the kinetic energy term. We focus in particular on values of the hopping terms which produce a flat band and show that, in the presence of contact and near neighbor repulsion, three phases exist: Mott insulator (MI), charge density wave (CDW), and the topological Haldane insulating (HI) phase which displays edge states and particle imbalance between the two ends of the system. We find that, even though the entanglement spectrum in the Haldane phase is not doubly degenerate, it is in excellent agreement with the entanglement spectrum of the Affleck-Kennedy-Lieb-Tasaki (AKLT) state built in the Wannier basis associated with the flat band. This emphasizes that the absence of degeneracy in the entanglement spectrum is not necessarily a signature of a non-topological phase, but rather that the (hidden) protecting symmetry involves non-local states. Finally, we also show that the HI phase is stable against small departure from flatness of the band but is destroyed for larger ones.Comment: 10 pages, 16 figure

Phase separation in the bosonic Hubbard model with ring exchange

We show that soft core bosons in two dimensions with a ring exchange term exhibit a tendency for phase separation. This observation suggests that the thermodynamic stability of normal bose liquid phases driven by ring exchange should be carefully examined.Comment: 4 pages, 6 figure

Solitons in Trapped Bose-Einstein condensates in one-dimensional optical lattices

We use Quantum Monte Carlo simulations to show the presence and study the properties of solitons in the one dimensional soft-core bosonic Hubbard model with near neighbor interaction in traps. We show that when the half-filled Charge Density Wave (CDW) phase is doped, solitons are produced and quasi long range order established. We discuss the implications of these results for the presence and robustness of this solitonic phase in Bose-Einstein Condensates (BEC) on one dimensional optical lattices in traps and study the associated excitation spectrum. The density profile exhibits the coexistence of Mott insulator, CDW, and superfluid regions.Comment: 5 pages, Latex with figure

Density of States and Magnetic Correlations at a Metal-Mott Insulator Interface

The possibility of novel behavior at interfaces between strongly and weakly correlated materials has come under increased study recently. In this paper, we use determinant Quantum Monte Carlo to determine the inter-penetration of metallic and Mott insulator physics across an interface in the two dimensional Hubbard Hamiltonian. We quantify the behavior of the density of states at the Fermi level and the short and long range antiferromagnetism as functions of the distance from the interface and with different interaction strength, temperature and hopping across the interface. Induced metallic behavior into the insulator is evident over several lattice spacings, whereas antiferromagnetic correlations remain small on the metallic side. At large interface hopping, singlets form between the two boundary layers, shielding the two systems from each other.Comment: 7 pages, 6 figure